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ABSTRACT 

Mrk 938 is a luminous infrared galaxy in the local Universe believed to be the remnant of a 
galaxy merger It shows a Seyfert 2 nucleus and intense star formation according to optical 
spectroscopic observations. We have studied this galaxy using new Herschel far-IR imaging 
data in addition to archival X-ray, UV, optical, near-IR and mid-IR data. Mid- and far-IR 
data are crucial to characterise the starburst contribution, allowing us to shed new light on its 
nature and to study the coexistence of AGN and starburst activity in the local Universe. The 
decomposition of the mid-IR Spitzer spectrum shows that the AGN bolometric contribution to 
the mid-IR and total infrared luminosity is small (Lboi(^GA^)/LiR ~0.02), which agrees with 
previous estimations. We have characterised the physical nature of its strong infrared emission 
and constrained it to a relatively compact emitting region of ^ 2 kpc. It is in this obscured 
region where most of the current star formation activity is taking place as expected for LIRGs. 
We have used Herschel imaging data for the first time to constrain the cold dust emission with 
unprecedented accuracy. We have fitted the integrated far-IR spectral energy distribution and 
derived the properties of the dust, obtaining a dust mass of 3 x 10^ M©. The far-IR is dominated 
by emission at 35 K, consistent with dust heated by the on-going star formation activity. 

Key words: Galaxies: evolution — Galaxies: nuclei — Galaxies: Seyfert — Galaxies: struc- 
ture — Infrared: galaxies — Galaxies: individual: Mrk938. 



1 INTRODUCTION 

Observations over the past decade have revealed that supermassive 
black holes (SMBH s) likely reside at the cent res of all galaxies with 
spheroids (see .Kormendv & Richston3l 19951 , for a review) and that 
the properties of these black holes and their host galaxies are tightly 

* Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 
ipation from NASA. 

** E-mail: pilar.esquej@cab.inta-csic.es 

* Augusto Gonzalez Linares Senior Research Fellow 



correlated (e.g. lMarconi & Hunlll2003l : lFerrarese & Forj|2005l) . In 
addition, the similarity in the anti-hierarchical evolution of both 
star formation rate and active galactic nucle i (AGN) activity at z<2 
jLa Franca et alj|2005l : lAmouts et al.ll2007h reveals that the assem- 
bly and evolution of galaxies and the SMBHs at their centres are 
intimately connected. However, the nature of such an intrinsic con- 
nection is still enigmatic. Since it is believed that the main phase 
of SMBH growth through major mergers occu rs in environments 
heavily obscured by dust (iHopkins et ai] |2006h . the infrared (IR) 
domain is the path to follow to unveil key signatures that remain 
lurking at other wavelengths. 
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The class of objects known as luminous infrared galaxies 
(LIRGs; LiR = 10" - 10'^ L^, IR in the range 8-1000 ^.m) 
contributes a significant fraction (^50%) of the cosmic infrared 
background. In addition, they are major contributors to the star- 
formation activity at z~ 1 felbaz et"al]|2002l : iLe Floc'h et^llOOSl : 
iPerez-Gonzalez et al. iHooir LIRGs are powered by star forma- 
tion and/or AGN activity, and a large fraction of them, especially 
at the high IR luminosity end, are classified as interacting galax- 
ies and mergers (Sanders & Mirabel 1996). It has been suggested 
that the mechanisms that trigger mass accretion onto the SMBH 
through major merger proce sses are episodes of nuclear star for- 
mation jHopkins etai]|2006h . Therefore, it is clear that studying 
LIRGs can shed new light on our understanding of the coeval evo- 
lution of host galaxy and black hole. 

Mrk 938 (also known as NGC 34, NGC 17, and 
IRAS F00085-1223, among others) is a local LIRG belonging 
to the IRAS Revised Bright Galaxy Sample ( Sanders et al. 20031) 
- RA=00hllm06.65s, dec=-12d06m26.7s. It is located at a dis- 
tance of 86.4 Mpc using the recession velocity relative to the Loca l 
Group czlg=5961 km s"' measured bv lSchweizer & Seitzed ^20071) 
and assuming a ACDM cosmology with (JIm. ^Ia) = (0.3, 0.7) 
and Hq = 70 km s"' Mpc"'. The measured IR luminosity of 
this galax y is Lir = 3 . 4 x 10 " Lq using the IRAS fluxes re- 
ported by ISanders et al. I ( I2OO3I ) and the adopted distance. Al- 
though this source has been classified in the literature as both 
Seyfert 2 and starburst (SB) galaxy from optical spectroscopy (see 
ISchweizer & Seltzer 2007. for a detailed discussion on this issue), 
its X-ray emission con firms the presence of an obscured AGN (see 
iGuainazzi et al. J| 2005L and also S ection 13. 11 1. We note that in the 
recent work of lYuan et alj ( l2010h the nuclear activity of Mrk 938 
has been definitively classified as Seyfert 2 based on optical spec- 
troscopy. 

The optical morphology of Mrk 938 is dominated by a bright 
red nucleus, a blue exponential disk, and a rich system of young 
massive star clusters. The faint optical surface brightness emis- 
sion reveals a pair of tidal tails to the northeast and southwest of 
the system, with the former being the brightest and longest one 
(see Fig. [H courtesy of F . Schweizer). All these signatures led 
ISchweizer & Seitze to propose that Mrk 938 is the remnant 

of a gas-rich merger of two unequal mass galaxies. The nuclei of the 
two galaxies have probably coalesced, creating a concentrated and 
obscured starburst and an AGN, both driving a strong gaseous out- 
flow. The starburst appears to be confi ned to a highly obscured cen - 
tral region of less than 1 kpc in radius jSchweizer & Seitzer*'2007'). 
This is indeed confirmed by ground-based high-angular resolution 
observations in the mid-IR. These revealed that the m id-IR emis- 
sion i n this system originates in the central ~ 800 pc jMiles et"al] 
Il996h . but it actually d ominates in a region of ~ 400 pc (see dis- 
cussion in Section |4.U . lFemandez et al ([20lrf) detected H i emis- 
sion from both tidal tails and from nearby galaxies, suggesting that 
Mrk 938 is part of a gas-rich group of galaxies. They also found a 
hint of emission between Mrk 938 and NGC 35, the largest com- 
panion lying at a projected distance of 131 kpc, and suggested that 
they might have recently interacted. The radio-continuum emis- 
sion detected in the nuclear region is extended, indicating that it 
is mostly due to the highly concentrated starburst. 

Here we present for the first time a detailed analysis of novel 
far-IR (70 - 500 yum) imaging observa tions of Mrk 938 obtained 
with the Herschel Infrared Observatory jPilbratt et al I20I0 ). These 
have been combined with our own analysis of X-ray, UV, and IR 
data already presented in the literature. The main goal of this work 
is to understand the AGN/SB connection in Mrk 938, as well as the 




30 kpc 



Figure 1. Ground-based 2.5' X 3' B band image of Mrk 938 provided by 
F. Schweizer. Note the bright nucleus and spiral stmcture in addition to the 
two tidal tails emerging from the system in opposite directions. 

origin of the IR emission of this galaxy. This is the second paper in 
a series presenting Herschel imaging observations of a sample of 
nearby Seyfert galaxies to investigate the co-existence of AGN and 
starburst activity, as well as the nature of the dusty torus. A study 
of the dust properties of NGC 308 1 has recently been presented in 
iRamos Almeida et alj ( 1201 Ibl) . 

This paper is structured as follows: Sect. 2 describes the data 
analysis procedure of the different IR and X-ray observations. In 
Sect. 3 we present evidence of AGN activity, including the decom- 
position of the mid-IR spectrum into AGN and SB contributions. 
The properties of the SB component have been derived in Sect. 4. 
We used the wealth of archival data to fit the SED of the galaxy 
and to derive the dust properties, presented in Sect. 5. Finally, our 
findings and conclusions are summarised in Sect. 6. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Herschel 

We performed far-IR imaging observations of Mrk 938 with 
the Photodetector Array Camera and Spectrometer (PACS, 
IPoglitsch etallbOld) and the Spect ral and Photometric Imaging 
REceiver (S PIRE, iGriffin et alj|20iol) onboard the Herschel Space 
Observatory dPilbratt et al1l2010l) . covering the spectral range 70 - 
500 (im in 6 bands. The data are part of the guaranteed time pro- 
posal "Herschel imaging photometry of nearby Seyfert galaxies: 
testing the coexistence of AGN and starburst activity and the na- 
ture of the dusty torus" (PI: M. Sanchez-Portal). 

The PACS observations were carried out using the "mini-map" 
mode, consisting of two concatenated scan line maps, at 70° and 
1 10° (in array coordinates), at a speed of 20"/sec, each one with 10 
lines of 3' length and cross-scan step of 4". This mode produces 
a highly homogeneous exposure map within the central 1' area. 
The set of maps were duplicated to observe through both the 
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70 \im ("blue") and 100 \im ("green") filters. Therefore the galaxy 
was observed twice through the 160 [im ("red") filter. The PACS 
beams at 70, 100, and 160 [im are 5.6", 6.8", and 11.3" FWHM, 
respectively. With the SPIRE photometer, the three available bands 
were observed simultaneously using the "small map" mode, with 
two 1x1 nearly orthogonal scan lines, at a scan speed of 30"/sec. 
The scan line length is 11.3' and the area for scientific use is around 
5'x5'. The FWHM beam sizes at 250, 350, and 500 \im are 18.1", 
25.2", and 36.9" respectively. 

The data reduction was carried out with the Herschel In- 
teractive Processing Environment (HIPE) v6.0.1951. For the 
PACS instrument, the extended source version of the standard 
PhotProject reduction script was deemed as adequate, due to the 
small angular size of the galaxy. This procedure implements a high- 
pass filtering algorithm to remove the 1/f noise of the bolometer 
signals. The sources are masked in order to prevent the high-pass 
filter to remove flux from extended stru ctures. The FM v5 pho- 
tometer response calibration files (Mtill er et alj|201lh were used. 
Level 2 maps produced by the reduction script for the two scan di- 
rections were merged using the mosaic task. For SPIRE, the stan- 
dard small map script with the 'naive' scan mapper task was ap- 
plied, using the calibration database v6.1. Colour corrections (for 
PACS, see Poglitsch et al. 2010; for SPIRE, please refer to the Ob- 
sen'er's Man ual) are small for blackbodies at the expe cted tem- 
peratures (e.g. lPerez Garcia & Rodriguez Espinosall200lh and have 
been neglected. More details on the processing of Herschel data are 
given in Sanchez-Portal et al. (in preparation). The far-IR maps of 
Mrk 938 are shown in Fig. [2] 

Aperture photometric data in the different bands were ex- 
tracted using IRAlQ and the measurements for an aperture radius 
of 70 "are shown in Table[T] The selected radius is a good compro- 
mise between including all IR emission and minimising the back- 
ground contribution. Since the PACS encircled energy radius of FM 
v5 is normalised to an aperture of 60", no aperture correction for 
the data was required. The morphology of Mrk 938, which will 
be further discussed in Sect. 14. 11 is compatible with an unresolved 
object. Using our far-IR and L 24um photometry and the parametri- 
sation of IPale & Heloi] ilOO'jt we derived a total IR luminosity o f 
LiR = 3.3 X 10" Lq, in agreement with that o f ISanders et aLH2003h . 

2.2 Spitzer 

2.2.1 Imaging Observations 

Mrk 938 was o bserved wit h the Spitzer i nstruments IRAC 
jFazio et al ."'2004') and MIPS (Rieke et al.'2004'). We retrieved the 
basic calibrated data (BCD) from the Spitzer archive (IRAC - 
P3269, PI: J.F Gallimore; MIPS - P3672, PI: J.M. Mazzarella). 
The BCD processing includes corrections for the instrumental re- 
sponse (e.g. pixel response linearisation), flagging of cosmic rays 
and saturated pixels, dark and flat fielding corrections, and flux cali- 
bration based on standard stars (see the IRAC and MIPS instrument 
handbooks for details). We combined the BCD images into mo- 
saics using the MOsaicker and Point source Extractor (MOPEX) 
software provided by the Spitzer Science Center (SSC) using the 
standard parameters. IRAC images of Mrk 938 at 3.6 and 8 (xm are 
shown in Fig. [3] The IRAC 8 |J.m image was affected by the banding 

' Image Reduction and Analysis Facility (IRAF) software is distributed by 
the National Optical Astronomy Observatories, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 



Table 1. Integrated Infrared SED of Mrk 938. 
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efi'ect (see IRAC Handbook). Similar to lOaUimore et~aLl ( |2010|) we 
corrected this artefact by fitting a polynomial to the affected rows. 

We computed integrated photometry of the source in the MIPS 
24 [xm image, and obtaine d from the literature photometric data 
points for the IRAC bands (Gallimore et al.l20IOl) . The source pho- 
tometry is presented in Table[T] 

2.2.2 Spectroscopic Observations 

We retrieved SpitzerflRS jHouck et alj|2004h spectroscopic obser- 
vations of Mrk 938 from the Spitzer archive (P3269, PI: J.F. Gal- 
limore; P30291, PI: G. Fazio). The observations were taken in low- 
resolution {R ~60-120) with the Short-Low (SL) and Long-Low 
(LL) modules and in high-resolution {R ~ 600) mode with the 
Short-High (SH) and Long-High (LH) modules. The low resolution 
data were observed in mapping mode while the high resolution data 
were obtained in staring mode. 

For the staring data we started with the BCD data. Bad 
pixels were corrected using the IDL package IRSCLEAlj^ Then 
we subtracted the sky emission and extracted the spectra using 
the standard programs included in the Spitzer IRS Custom Ex- 
traction (SPICE) package provided by the SSC and the point 
so urce calibration. The re sulting spectrum was already presented 
in iTommasin et al] ( l20I(lt) . For the SL+LL ma pping observation s 
we constructed the data cubes using CUBISM ( Smith et"al]|2007l) . 
We extracted the nuclear spectrum of Mrk 938 observed in SL-l-LL 
mapping mode using a 13'.'4xl3'.'4 aperture. Since the mapping 
mode data cubes are calibrated as extended sources we applied 
a wavelength dependent aperture correction to obtain a point-like 

^ The IRSCLEAN package is available at 

|http://irsa.ipac.caltech.edu/data/SPITZER/docs/dataanalysistools/tools/irsclean/| 
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Figure 2. Herschel images of 3' X 3' centred on Mrk 938. North/East is up/left respectively. Top: PACS 70, 100 and 160 i^m; bottom: SPIRE 250, 350 and 
500 \aa\ from left to right respectively. The PACS images show the shape of the PSF. In particular, a plume-like structure lying approximately in the NE-SW 
direction has been deemed as a PSF artifact and it should not be confused with the tidal tail observed in the IRAC images (Fig. [3). 



spectrum to represent the nuclear emissi on of this galaxy. For 
more d etai ls on the data analysis please see lPereira-Santaella et al.l 
( l2010ij|bl) : lAlonso-Herrero et al.l j2012h . 

We evaluated the strength of the 9.7 \im silicate fea- 
ture, 5 Si, by calculating the ratio of the observed flux den- 
sit y to the continu um flu x densit y at 9.7 ii.m using the method 
in IPereira- Santaell a et alj j2010b[) . Applying the definition of 
ISpoon et"al.l ( I2007D we get 5si = -1-1, where the minus 
sign indicates absorption. This value is considerably higher, 
in absolute terms, than th e average of the LIRG sample in 
lAlonso-Herrero et al.l ( 1201 2h of (5si> = -0.6. We also measured 
the equivalent width (EW) of th e 6.2 |i.m PAH featur e, and obtained 
EW(PAH 6.2 ^m)=0.4±0.1 um. lSpoon et"al] ( |2007|) provided a di- 
agnostic diagram to give a general classification of IR galaxies 
based on 5si and EW(PAH 6.2 \im). According to that, Mrk 938 
is a moderately obscured composite nucleus. 

2.3 XMM-Newton 

Mrk 938 was observed by XMM-Newton on 2002, December 22 
(rev 556, obsid 0150480501, PI: R. Maiolino) and then anal- 
ysed with the XM M-Newton Science Analysis Software (SAS; 
iGabriel et Zl |2004|) . All EPIC observations were taken in Full 
Frame mode and the medium filter in place. These data were al- 
ready presented in Guainazzi et al. (2005), but we reduced them 
with the latest software to take full advantage of the most recent 
calibration files. 

To extract the X-ray spectra, high background periods pro- 
duced by intense soft proton fluxes were rejected by defining good 



time intervals within the complete observation window. After this 
process, the net exposure times were about 13, 16 and 16 ks for pn, 
MOS 1 and M0S2 respectively. The light curves of the sources were 
inspected and no variability within the X-ray observations was ob- 
served. Source photons were extracted from a circular region of 30 
arcsec radius centred on the object position. A circular source-free 
region on the same chip and radius of 60 arcsec was used to char- 
acterise the background. The target does not suffer from pile-up, so 
the pn spectrum was extracted with patterns to 4, and MOS spec- 
tra from to 12. Ancillary files and response matrices were gen- 
erated by the arfgen and rmfgen SAS tasks respectively in order 
to convert the counts to physical units during the spectral analy- 
sis. All EPIC spectra were binned to oversample the instrumental 
resolution and to have no less than 20 counts in each background- 
subtracted spectral channel. 

Spectral fitting was performed with XSPEC 12, where all ap- 
plied models include Galactic foreground abs orption (A^h=2.61 x 
10^° cm^^) inferred from H i observations dOickev & Lockmiuil 
1 19901) . Given that our source has been found to be highly obscured, 
we started fitting the data with a model consisting of a powerlaw 
plus an absorbed powerlaw, both with the same photon index. We 
got a reasonable fit at high energies, but a prominent soft excess 
with respect to the powerlaw appeared below ~ 1 keV. Thus we 
added a thermal component, MEKAL in XSPEC. We obtained a tem- 
perature ^:r=0.67keV typical of a starburst emission, and a high 
intrinsic absorption of N'^'=15*3,\ x 10^" cm"^, consistent with high 
obscuration. A considerably flat powerlaw was obtained, rx=1.43. 
A summary of the parameters of this model, which is considered 
our best-fit (^f^/dof=50/64), is presented in Table (2] Our derived 
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Table 2. Parameters derived from the fitting of tlie best-fit model, 
zwabsx(mekal+zpo+zwabsxzpo) in XSPEC nomenclature, to the XMM- 
Newton data. Galactic foreground absorption has been used for the first 
model component. Luminosities have been corrected for absorption. Errors 
are 90% confidence. 
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Figure 3. SpitzerflRAC images at 3.6 (top) and 8 ]im (bottom) of Mrk 938 
covering a field of view of 3' X 3'. The orientation is north up, east to the 
left. Images are shown in a square root scale. Note that the northeast tidal 
tail and knots of star formation are clearly detected in the 3.6 |j.m image, 
which traces photospheric emission. In the 8 i^m image, which probes dust 
continuun and PAH emission, one can see faint emission in that tail likely 
to be associated with the knots of star formation also seen in the optical (see 
Fig.EJ. 



■0^ 


11 

. . 1 













"l 2 ' ' 5 

channel energy (keV) 



Figure 4. X-ray spectrum of the EPlC-pn data. The upper panel shows data 
and best-fit model (solid line). Residuals in terms of sigmas are plotted in 
the bottom panel. 

2.4 Swift 

The position of Mrk 938 has been obs erved with the Burst Alert 
Telescope (BAT iBarthelmv et ahlllOOS i). but the source is not de- 
tected even in the 70th month survey. We have used the stan- 
dard BAT tool batcelldetectl^] to analyse the observations and 
derive the rate (normalised by the number of working detectors) 
at the source position in the eight standard channels from cor- 
responding background subtracted images. In the fourth channel, 
35-50 keV, the rate becomes negative, therefore we have com- 
bined the first three channels to get a rate in the 14-35 keV en- 
ergy range. Using WebPIMM^ for the count rate to flux conver- 
sion we get a flux upper limit of 2.5xl0"'2erg s"' cm"^ (equiv- 
alent to a luminosity of 2.2x10"*- erg s"'). If we extrapolate our 
best-fit model derived for the XMM-Newton data, we get Fi4_35ijev 
= 1.6+'^qI X 10-'- erg s"' cm-2(L,4-35keV =l-4!o;5 xlO"*- erg s"'), that 
is compatible with a BAT non-detection. We cannot further con- 
strain the X-ray model because the extrapolated fluxes for all rea- 
sonable models lie below the BAT upper limit. 



X-ray luminosity is similar to that in iGuainazzi et"d] j2005h . Fig. 
|4] shows the X-ray data (only EPIC-pn data is shown for simplic- 
ity) and model-t o-data residuals. We also tried adding a reflection 
model (PEXRAV, iMagdziarz & Zdziarskilll995h . for possible X-ray 
scattering from an optically thick surface, but did not get any sta- 
tistical improvement in the fit. 



3 AGN ACTIVITY 

In the local universe, a significant fraction of LIRGs hosts an 
AGN, and this fraction increases with the IR luminosity of the 
system tVeilleux et al., ,1995i) . Using optical emission-line ratios 



3 http://heasarc.gsfc.nasa.gov/ftools/caldb/help/batcelldetect.htmll 
^ http://heasarc.nasa.gov/Tools/w3pimms.htmll 
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and a new scheme to classify galaxies, lYuan et ^ ( I2OI0I) re- 
cently reported that approximately half of local LIRGs are clas- 
sified as Seyfert, LINER or composite (SB/AGN). The AON de- 
tection rate in LIRGs can be as high as ~ 60% when combin- 
ing optical and IR indicato rs (i.e. emission line ratio d iagnostics 
and spectral decomposition. [Alonso-Herrero et aLll2012h . and thus 
similar to that of l ocal ultraluminous i nfrared galaxies (ULIRGs, 
LiR = 10'- - 10''' Lm lNardini et al.l20I O'). However, in clear contrast 
with ULIRGs, in local LIRGs the AGN bolometric contribution to 
the IR luminosity of these sys tems is very small, with an average 
of <Lboi(AGN)/LiR>=0.05!jJ°^ jAlonso-Herrero et al.ll2012l ). There- 
fore the high IR luminosities in mos t LIRGs are produced in intense 
dusty starbursts. lPetric et al.l 1 I2OI H) estimated that AGN contribute 
~12% to the IR luminosity, wh i ch is c ompatible with the upper end 
found in lAlonso-Herrero et al. ||20l2). 

3.1 X-ray evidence 

The analysis of the X-ray emission of Mrk 938 shows evidence 
of the presence of an obscured AGN. The active nucleus domi- 
nates the hard X-ray contribution as can be seen from the highly 
absorbed powerlaw component used to fit the spectrum. The mea- 
sured absorption-corrected hard X-ray luminosity of 1-4^552 x 



|4.3l l. The high intrinsic absorbing column returned by the spectral 
fitting supports the conclusion of the highly obscured nucleus. 

A powerful diagnostic to estimate the nuclear obscuration is 
the 2-10 key to [Om] flux ratio, 7 (e. g. iMaiolino et ^Il998l : 
ISassani et al.lll999l : lGuainazzi et al.ll2005l) . According to previous 
studies the threshold lies in r=0.1, where objects with T ^0.1 
are Compton-thick while those with T ^1 are unobscured or 



C ompton-thin. Using F[o mi = 1.1x10 '''erg s ' 



cm as presented 



m IPolletta etaLl ( ll99d) we obtain T=0.01 for our source, therefore 
Mrk 938 would appear as Compton-thick. However, this galaxy has 
a very strong starburst contribution in the optical (Goncalves et al.l 
1 199^ ). therefore the AGN is not responsible for a considerable frac- 
tion of the [O III] flux and the T ratio as a measure of obscuration is 
uncertain. Another indicator of the nuclear obscuration is the emis- 
sion fluorescence Fe K„ line. Objects are considered as Compton 
thick when EW(Fe K„)>1000eV (Matt et al. 1996). In the case of 
Mrk 938, due to low statistics at the hardest energies we can only 
constrain the presence of t he Fe Kg line by extra cting an upper 
limit. Using Cash statistics, iGuainazzi et alj ( |2005|) fitted the un- 
binned spectrum in the 5.25-7.25 keV range. They found an equiv- 
alent width upper limit of 322 eV when fitting a Gaussian with en- 
ergy fixed at 6.4 keV. Therefore, the width of the Fe K„ line places 
the source in the Compton-thin region. 



3.2 Mid-infrared evidence 

In terms of continuum emission, the mid-IR spectral range con- 
tains a wealth of information to disentangle the AGN from the 
SB contribution in galaxies. This comes from both broad fea- 
tures and fine structure emission lines. In particular, high metal- 
licity starburst galaxies show very similar mi d-IR spectra (see e.g., 
iBrandl et alj2006l : lBemard-Salas et al.l2009b . which are rather dif- 
ferent from the continuum and line emission observed in pure 
AGN. Because of this, spectral decomposition of Spitzer IRS spec- 
tra is a powerful tool to estimate the fractional AGN and SB con- 
tributions at d ifferent mid-IR wavel engths in both loc al (see e.g., 
iNardini et al1l20 08. 2010; Alonso- Herrero et al.ll2012h and high-z 
(e.g. Menendez-Delmestre et al.ll2009h galaxies. 



Table 3. AGN fractional contributions within the IRS slit obtained from the 
AGN+SB decomposition. 



C^mfAGN] C'f^^lAGN] C'f^^[AGN] C^f^^[AGN] 



0.21 



0.10 



0.05 



0.02 



In the case of local LIRGs lAlonso-Herrero et al.l j2012h 
demonstrated that the mid-IR spectral decomposition can identify 
even subtle AGN emission that may otherwise be under the de- 
tection threshold with other mid-IR diagnostics. Here we used a 
spectral decomposition method very similar to theirs to estimate 
the AGN contribution with respect to the total mid-IR emission 
of Mrk 938. For the SB component this me thod uses the aver- 
age starburst spectrum o f' BrandletaL * (200^), and the templates 
of purely star forming (U)LIRGs oflRieke et al. (200.2) covering 
the IR luminosity range 10.5 < log(LiR/Lo) ^ 12. The AGN 
continuum emission is represented by that prod uced by an in- 
terpolated version of the CLUMPY torus models of iNenkova et al.l 
('2008a, b). We restricted the large database of the CLUMPY torus 
models to only use two AGN templates. These were derived by 
iRamos Almeida et alj ( 1201 lal) as the best fitted models to the aver- 
age S EDs of nearby type 1 and type 2 Seyferts using a Bayesian 
code dAsensio Ramos & Ramos Alme ida 2009). We then looked 
for the best combination of SB -1- AGN templates to fit the IRS 
spectrum that minimised x^- The best fit to the IRS spectrum of 
Mrk 938 was achieved with the type 2 AGN template, as expected, 
plus the 10'^ Lq starburst template. This SB template is required to 
fit the relatively deep 9.7 [im silicate feature observed in this galaxy. 
We note that the IR luminosity of the best fit template is higher than 
that of our galaxy. This might indicate that the IR emission comes 
from a relatively compact region (Sect. l4Tt . as is the case of many 
LIRGs. 

In Fig. |5] we show the IRS spectrum of Mrk 938 plus the 
AGN-l-SB best fit model, in addition to the scaled model compo- 
nents. As can be seen from the figure, the AGN contribution to the 
mid-IR emission is small but detectable. The AGN fractional con- 
tributions to the continuum emission within the Spitzer/IRS slit are 
listed in Table [3] Although our statistical uncertainties on the de- 
rived parameters are small (±5%) the errors associated to the tem- 
plates can be larger, especially in th e SB comp onent. Our results 
are in good agreement with those of lOeo et al.l ([2009^. Clearly, as 
found in local LIRGs and many local Seyfert galaxies, the AGN 
cont ribution at wavelengths longer than 40 \im is gre a tly reduced 



cont ribution at wavelengths longer tnan 4U \xm is gre a tly reducea 
("see I Alonso-Herrero et aLll2012l ; lMullanev et al.ll201 ll) . IVega et al.l 
( I2OO8I) obtained a contribution of 17% of the AGN to the 3- 
30 vim luminosity, which is at the upper end of th e AGN contribu- 
tion th at we find at 20nm including uncertainties. iTommasin et al.l 
( l2010h find an AGN contribution to the I9[un flux of Mrk 938 of 
50(±I3)%, which is considerably higher than our estimation. This 
can be explained because they used upper limits on the undetected 
line widths for their estimation. 

Additionally, we can use mid-IR high excitation emission 
lines to investigate the AGN emission, in particular, the [Ne v] 
lines at 14.32 yim a nd 24.32 yim and the [O iv] /125.89 (im line 
dGenzel et ai]|l998l) . lTbmmasin et al.l ( l2010h did not d etect any of 
these emission lines in Mrk 938 (see also iGallimor e et al. 201(3; 
IPereira-Santaella et alj|2010ah . When these high excitation emis- 
sion lines are ratioed with a line mostly produced by star formation 
(e.g., the [Neil] /112.81 [im), it is possible to obtain an approxi- 
mate contribution -if any- of the AGN to the mid-IR emission of a 
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Figure 5. AGN+SB best fit of the observed SL+LL IRS spectrum of 
Mrk 938 (black line). The dashed-dotted green line is the fitted scaled 
Seyfert 2 torus model template (see text for details), the purple dashed line 
is the fitted scaled SB template, and the solid red line is the sum of the fitted 
SB and AGN components. 



galaxy. The upper limit to the [Oi v ] /125.89 vim/[Nen] /tl2.81 |^m 
of ~ 0.01 jTommasin et alj|201(]l : IPereira-Santaella et aD l2010ah 
indicates that the mid-IR emission of this galaxy is dominated 
by star formation activity (see e.g., lAlonso- Herreroetal]|2012l : 
IPetric et aDl201 ll) . This is in good agreement with the low AGN 
fractional contribution to the mid-IR emission inferred from the 
spectral decomposition. 

Since the IRS slit encloses most of the 24 (i.m emission (from 
the comparison with the MIPS photometry, see Table [T), the de- 
rived AGN fraction contrib ution to the 24 \xm emission is 0.05. 
lAlonso-Herrero et al ] ( l2012h also found that in local LIRGs this 
fractional AGN contribution is related to the AGN bolometric con- 
tribution to the IR luminosity of the system. Using their relation 
we estimated that in Mrk 938 Lt,„,(AGN)/LiR = 0.02+j;o2 xhis 
is in good agreement with th e AGN bolometric contribution es- 
timated by IVega et al] ( 1200 8l) of 4% from their modelling of the 
full SED of this galaxy. We derived then an AGN bolometric lumi- 
nosity of 2.6^^3 xlO'*' erg s"'. Alternatively, the AGN bolometric 
luminosity can be estimated from the hard X-ray 2 - lOkeV lu- 
minosity. Using the bolometric correction of lMarconi et al.l ( l2004h 
and our derived X-ray luminosity from the model fit L2-iokev = 
1.4 X lO'^^erg s"' (Table |2j, we get a bolometric luminosity of 
Li,oi(AGN)=1.6^Qg xlO*^^ erg s"', compatible with our previous es- 
timation derived from the IRS spectral decomposition. 



4 STARBURST ACTIVITY 

It is now becoming clear that the size of star-forming regions and 
the resulting IR surface brightness (also referred as compactness) 
of IR bright galaxies may help us understanding how star formation 
is produced in galaxies. That implies that we may be able to differ- 
entiate between extended steady star formation in disks and com- 
pact star formation activity driven by mergers. These two processes 
produce different overall IR spectral energy distributions (SEDs) 
of LIRGs and ULIRGs both locally and at high-redshift (see e.g., 
iRuiopakarn et al.ll2011al lbl and references therein), where there is 
no clear dependence of the shape of SED with Lir. 

Similar to the fundamental planes of spiral and elliptical 
galaxies, there are three parameters that define a fundamental plane 



of IR-selected galaxies. ( IChanialetai]|2007l) . These are: the IR 
luminosity, the dust temperature and the size of the IR-emitting 
regions. According to these authors, changes in the IR compact- 
ness describe a smooth sequence going from quiescent to starburst 
galaxies, where there is a simultaneous increase of gas surface den- 
sity, effective dust temperature, and SFR surface density. In the fol- 
lowing sections we explore in detail some of these properties in 
Mrk 938. 



4.1 IR Emitting Region 

The detailed study of the optical properties of Mrk 938 (mor- 
phology, star cluster ages, and dynamics) led Schweizer & Seitze^ 
( l2007t) to claim that this object is the result of a merger of two gas- 
rich galaxies of unequal masses (l/3< m/M <2/3). At present, the 
current starburst activity of Mrk 938 appears to be confined to a 
highly obscured central region of less than 1 kpc in size. 

Fig. |3]presents the SpitzerfiKAC images at 3.6 and 8 (^m. The 
3.6 (im image, which probes the photospheric emission, shows the 
main body of the galaxy plus some extended low surface brightness 
emission. The northeast tidal tail of this merger system is clearly 
detected in the form of diffuse emission and some knots. There also 
exists a southwest tail with a much fainter optical surface brightness 
that is hardly seen in the 3.6 vim image. The knots in the northeast 
tail, also detected in the IRAC 8 |xm image, appear to be coinci- 
dent with the blue knots lyi ng at a projected distance of ~ 50" 
(~ 21 kpc) from the nucleus jSchweizer & Seitzeill2007h . Knots in 
both tidal tails are associated with H n regions (A. Perez-Garcfa, 
private communication). 

The IRAC 8 vim emission, which probes dust continuum emis- 
sion and PAH emission, is dominated by the presence of a bright 
central region. This region appears slightly resolved with a size 
(FWHM) of 2.3" (~ 950 pc). The size of the 8 ^m emitting re- 
gion inferred from the IRAC data is in good agreement with the 
distribution of the mid-IR emission detected from g round-based 
high angular resolution imaging as measured bv lMiles e t al. ( 199^. 
These authors also showed that there are two mid-IR sources in 
the central regions of Mrk 938 separated 1.2", although most of 
the mid-IR emission arises from a resolved source with an angular 
size of 1" (FWHM). This confirms that the starburst region in this 
galaxy is only relatively compact (~ 0. 5 - 1 kpc, FWHM), as is the 
case in a large fraction of local LIRGs jAlonso-Herrero et alj2006l : 
IRuiopakarn et al.ll2011ah . Also, the presence of extended mid-IR 
emission in the central regions of Mrk 938 is in agreement with the 
result that the AGN contribution to the mid-IR emission is low, as 
demonstrated in Sect.[3i2] We used PACS 70 \im and MIPS 24 ^m 
images, that have comparable angular resolution, and probed the 
current star formation. We can set an upper limit to the starburst 
region of ^ 5" ~ 2 kpc (FWHM), which is in good agreement with 
the measured size at shorter wavelengths. 

We analysed the Herschel brightness profiles of Mrk 938 and 
compared them with those of point-like sources. We found no ev- 
idence of extended emission in any of the PACS or SPIRE bands. 
As a sanity check, we performed aperture photometry on the PACS 
70 vim image using small radii, in this case 6 and 14 ". After apply- 
ing the corresponding aperture corrections we obtained values in 
agreement with the integrated flux at 70 nm, therefore confirming 
that we are not detecting extended emission. 
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4.2 Star Formation Rate 

Based on a varie ty of morpliological, photometric and spectro- 
scopic arguments, ISchweizer & SeitzeJ (120071) concluded that the 
star formation activity of Mrk 938 was higher and more spread out 
in the past. As a result of the merger process a galaxy-wide star- 
burst started about 600 million years ago. It is this post-starburst 
population that appears to dominate the emi ssion from the blue 
disk present in this system. iRiffel et al] ( 1200 St) used near-IR spec- 
troscopy to study the stellar populations in the central region of 
this galaxy. They found that the near-IR emission of this galaxy 
is dominated by an intermediate-age stellar population of ~ 1 Gyr, 
whereas \he young 30Myr) stellar population contributes to less 
than 25% of this emission. 

In the previous section we saw that most of the IR emission 
of Mrk 938 appears to be coming from a relatively compact re- 
gion 0.5-2 kpc in size. In very dusty galaxies, as is the case of 
Mrk 938, this IR emission is probing the current star formation rate 
(SFR). This is because the UV light emitted by the young stars is 
absorbed by the dust and re-emitted in the IR. Therefore, some IR 
monochromatic luminosity or, al ternatively, the total IR luminosity 
can be used to estimate the SFR dKennicutt et al J 20091) . The 24 [im 
monochromatic luminosity has been fo und to be a good indicator of 
the current SFR of dusty galaxies (e.g . lAlonso- Herrero et alj2006l ; 
ICalzetti et ar]|2007l : iRieke et alj|2009l) . It has the advantage that it 
traces dust heated by the most massive stars, he nce directly probing 
the current star formation activity. We used the lRieke et al.l ^20091) 
calibration 



SFR,R(Moyr-') = 7.8 x 10-'%4pm (io) 



x{7.76x 10-"L24™(Lo)l' 



which assumes a Kroupa IMljfl Using L24^m = 5.3 x 
10'° Lq (calculated from the 24 \im flux in Table [T]l, we estimated 
a SFRiR of 42 Mq yr"' to an accuracy within 0.2 dex (after remov- 
ing the AGN contribution at this wavelength). This value is actually 
the obscured SFR . This value lies between the values estimated by 
IVega et alj ( l2008h of 21 and 15 1 Mq yr"' based on SED fitting for 
star formation in the last 10 Myr and averaged over the entire burst 
respectively. Also, our derived SFR is similar to the 64 M q yr"' 
obtained from 1.4 GHz observations dFemandez et alj|201(]|) . As a 
sanity check we can test the above hypothesis that most of the on- 
going SFR in LIRGs is taking place in obscured environments. 

A source of uncertainty in the SFR calculation using IR indi- 
cators arises if a significant fraction of the stellar luminosity is not 
absorbed and then re-emitted in the IR, but escapes directly in the 
UV instead. We derived the SFRuv using the far-UV flux for this 
galaxy. This was found by correlation of the source position with 
the GALEX archive. Taking into account a Galactic colour excess 
E{B - V)=0.027 (quantity reported in NEE@) and applying the ex- 
tinction curve of iFitzpatricS j 1 9991) we derived the UV corrected 
flux (F;jy=0.62mJy). FoUowing lKennicuttl h998l) . assuming a flat 
continuum and scaling to a Kroupa IMF we have 



0.048 



(1) 



SFRl 



,(Moyr-'): 



(ergs ') 



(2) 



Using a Salpeter IMF would result in a SFR a pproximately 1.4 times 
higher than that from the Kroupa IMF (see e.g. Ken nicutt et al.l200^ . 
^ The NASA/IPAC Extragalactic Database (NED) is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration 



where Ly is the UV luminosity. For Mrk 938 we find that the 
SFRuv represents ~3% of the SFRir, as also found in lHowell et al.l 
( l2010l) . This value is in agreement with the medi an value found for 
a samp le of LIRGs in Buat et al. ( 2007) and Pe reira-Santaella et al.l 
( I2OIII) . 

Similarly, we can estimate the unobscured SFR from the 
observed (not corrected for extin ction) integrated Ha luminosit y 
(Lhi, = lO"*- erg s"') reported by Moustakas & Kennicutt (200f^. 
We used the SFR recipe of . Kennicutt et al. (2009) that accounts for 
the unobscured and obscured SFR in terms of the Ha and 24 yaa 
luminosities, respectively. We inferred that most (~ 90 - 95%) of 
the current star formation in Mrk 938 is obscured, and thus likely 
to be taking place in the compact starburst region revealed by the 
mid-IR observations described above. This is in agreement with the 
small value found for the SFRuv calculated above. 



4.3 X-rays from star formation 

For a starburst, the emission of radiation in the X-ray domain 
can be attributed to a combination of components like super- 
nova remnants, hot plasmas associated to star-forming regions and 
High Mass X-ray Binaries (HMXB) within the galaxy among oth- 
ers (jpabbiano 1989)- Several investigations have revealed correla- 
tion between the X-ray emission and total SFR in starbursts (e.g. 
iRanalli et alj|2003l: IPereira-Santaella et al.ll20I ih . The linear rela- 
tion of lPereira-Santaella et alj j20I ll) reads 



SFRiR+uv [Mq yr ') = constant x 10 ^Lx erg s 



(3) 



with a 0.25 dex scatter, where constant='iAi±Q3) and 3.9 (±0.4) 
for the soft 0.5-2 keV and hard 2-10 keV X-rays respectively. Ap- 
plying this equation with the SFRs derived in Sect. 14.21 we esti- 
mate the following luminosities due to star formation activity in 
Mrk 938: Lo.5-2keV = 1.2(±0.6) x 10*' erg s^'and Lj-iokeV = 
1.0(±0.6) X lO'" erg s"'. The errors take into account the uncer- 
tainties in the SFR value and the uncertainties in Eq.[3] 

We have compared the X-ray emission expected from star for- 
mation activity as calculated above with that of the fitted XMM- 
Newton data. Using the luminosities measured in Sect. 12.31 for the 
best-fit model to the X-ray data we get a consistent value for the 
soft component (0.5;^q2 x 10"" erg s"'). This is in agreement with 
the fact that most, if not all, of the X-ray emission at soft energies 
is produced by the starburs t, as is also the ca se in a large fraction 
of local LIRGs (Pereira-Sa ntaella et al.ll20lll) . Also, the model de- 
rived from the fit of the X-ray data is consistent with star formation. 
However, in the 2-10 keV band, the luminosity of the X-ray fitting 
is an order of magnitude higher than that expected from the SFR. 
This was a priori anticipated because the AGN is the major con- 
tributor to the harder energies. 



5 SPECTRAL ENERGY DISTRIBUTION AND DUST 
PROPERTIES 

Several attempts have been made in the past to model the IR 
and submillimeter spectral energy distributions (SED) of LIRGs 
and ULIRGs using a modified blackbody with one or sev- 
eral dust temperature or full radiative trans f er modelling (see 



Klaas et al.l I2OO1I : iDunne & EaQ 1200 iL IVega etal.l 1200^: 



[elements et al, II2OIOI) . Here we use for the first time Herschel pho- 
tometry and data from the literature to obtain better estimates of 
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the dust properties of Mrk 938. Table[T]presents our Herschel pho- 
tometry as well as IRAS, Spitzer and ISO data compiled from the 
literature to construct the far-IR SED. The errors quoted for the 
Herschel flux densities take into account the photometric calibra- 
tion uncertainties of ~ 10% in all the bands. The remaining errors 
listed in Table[T]are those given in the corresponding references. 

We parametrised the SED of Mrk 938 using the standard mod- 
ified blackbody, which depends on the dust temperature (rju,,) and 
the dust emissivity (fi). In the case of optically thin emission in the 
IR, this function can be approximated as: 



(4) 



We fixed P = 2 as it has been found to be a reasonable value 
for IR bright galaxies llDunne & Ealesll200lh . Then we used a x' 
minimization method to fit the dust temperature. The SED fit was 
done using all the flux densities for A ^ 60 fxm given in Table [T] 
with the normalization at 100 \im. For the IRAS and ISO flux den- 
sities we used a 10% uncertainty. 

We obtained a best-fit dust temperature of Tdust = 36.5 K. As 
can be seen from Fig. [6] (left panel), using a single temperature 
produces a good fit to the far-IR SED of Mrk 938. Acceptable fits 
(X^ < 1.5 x Xmiii^ were also obtained in the range of temperatures 
T^dust = 36^3 K (the fits for the limitting temperatures are shown in 
Fig. ID as dotted lines) . The best-fit temperature is similar to those 
found by Klaas et al] ( 1200 ll) for local (U)LIRGs using similar as- 
sumptions. We also tried a fit using /3 = 1.5, which has been found 
to fit better the SEDs of LIRGs and ULIRGs founne & Ealesl200ll : 
[elements et"ai]|2010h . The fit provided a higher dust temperature, 
as expected. However, this fit was considerably worse, having a 
value of more than three times higher than that obtained with 
13 = 2. 

We derived the dust mass (Mdust) using the following equation 
(adapted from .Hildebrand 1983 ): 



KyBy(TiuSf) 



(5) 



where fy is the observed flux density, is the luminosity distance, 
ByiTiiist) is the blackbody emission for the best-fit dust temperature, 
and Ky is the absorption coefficient. We evaluated this expression at 
250 [im using an absorption coefficient of <C250um = 4 .99 cm ^ g"' 
as interpolated from the dust model of Li & Draind ( 1200 ll) . We 
obtained A/dust = 3 x 10^ Mq, which is consistent with the val- 
ues derived for l ocal ULIRGs and ot her IR-bright galaxies (see 
iKlaas et alj200ll : lciements et al.l2O10l and references therein). For 
the range of acceptable dust temperat ures, the inferred u ncertainty 
in the calculated dust mass is ~ 20%. IVega et al.l ( f2008l) estimated 
an Mdust of 8.5 x 10' Mg. However, they used the dust-to-mass ra- 
tio as conversion factor between the gas and the dust mass, using 
a fixed value of 100 for all sources in their sample. This can intro- 
duce a large uncertainty in the dust mass. This is most p robably the 
reason for the difference between our value and that in lVega et al.l 
( l2008l) . 

We have also fitted the SED of LIRGs using a t wo-temperature 
model, similar to those used in other works (e.g. iDurme & EalesI 
[200 1; Clements et al. 2010), where we have taken into account all 
flux density values for A ^ 24 [im to obtain a better constraint of the 
properties of the warm dust. Using two modified blackbody func- 
tions and j8=2 we obtained T^^^ = 67 K and T^^^^ = 35 K for the 
warm and cold components, respectively. As expected, the temper- 



ature of the cold dust is very similar to that of the single black- 
body fit. We also note that the temperature of the warm dust is 
not tightly constrained, as T^^^^^ = 67 ± 5 K produces values of x^ 
similar to;^^^j^. Finally, the goodness of the fit of the two tempera- 
ture model is only ma rginally better than the single temperature fit. 
[elements et alj ( 1201 Ol) found higher dust temperatures for ULIRGs 
than those of lower luminosity systems, i.e. a median dust temper- 
ature of 42 K for a sample of ULIRG s, while a m edian of 35 K has 
been found for normal galaxies tDunne & Ealesi .2001). The latter 
is consistent with the relatively cold temperature found in Mrk 938. 
We also find that the mass ratio of the cold to the warm dust com- 
ponent is very high (M^t^sJ^dust ~ which clearly indicates that 
the cold component is dominant in this galaxy. The total dust mass 
is very similar to the value obtained with one dust temperature. 

In Fig.|6](right panel) we have plotted the entire SED includ- 
ing the torus model of hot dust and the best-fit model to the far-IR 
data using the combination of two blackbodies described above. 
The latter best-fit temperatures are typical of star forming regions 
in the galaxy disk, therefore confirming that the active nucleus is 
not responsible for most of the observed dust emission (as seen in 
e.g. lHatziminaoglou et al.l l201Cf). Dust in regions of star formation 
have peak temperatures between 40 and 70 K, while radiation of 
AGN origin has a characteristic peak temperature within the range 
120-170 K (e.g. IPerez Garcia & Rodriguez EsDinosall200lh . This 
is in agreement with what we found in Sect. |3.2| from the spectral 
decomposition, where we concluded that the dusty torus only con- 
tributes up to 2% of the observed IR emission . 

As a comparison, iRamos Almeida et al. I( l2011bl) studied the 
dust distribution of the Seyfert 2 galaxy NGC 3081 also includ- 
ing Herschel PACS and SPIRE data in their analysis. The nu- 
clear emission of this galaxy was reproduced by a clumpy torus 
model of warm dust heated by the AGN. The circumnuclear emis- 
sion was fitted by a cold dust component at 28 K heated by young 
stars in the galaxy disk, plus a colder component at 19 K in the 
outsk irts of the galaxy heat ed by the interstellar radiation. Simi- 
larly, IRadovich et alj ( Il999 ) decomposed the SED of the narrow 
line galaxy NGC 7582 and found an extranuclear emission domi- 
nated by cold dust at 30 K, plus an additional very cold component 
emitting at 17 K. We do not detect this very cold dust component 
in Mrk 938. This is because we cannot resolve any extended emis- 
sion in our data, and although this component is likely to exist its 
contribution to the total IR luminosity is negligible. 

We can finally estim ate the size of the IR emitting region. Fol- 
lowing iKlaas etalJ (l200l'). for an optically thin (transparent) black- 
body with p = 2, the far-IR emitting component is homogeneously 
distributed in a "minimum face-on disk" of radius r^ with the form: 



'•t(pc) 



Mdust /Mq 



(6) 



i^41.7:/rrioo,/m ) 

Assuming an opacity at 100 [xm Tioo^m ~ 0.5 - 1, we obtain a 
minimum radius for the IR emitting region of ~ 500 - 800 pc, in 
good agreement with our measurements from the mid-IR images in 
Sect. dm. 



6 SUMMARY AND CONCLUSIONS 

LIRGs have emerged as an important cosmological class given that 
they may play a central role in our understanding of the general 
evolution of galaxies and black holes. We have studied Mrk 938, 
a LIRG in the local Universe classified as a Seyfert 2 through 



10 p. Esquej et al. 




Wavelength [/^m] Wavelength [/tim] 



Figure 6. SEDs of Mrk 938. The left panel shows the far-IR SED and the best-fit (solid line) with a modified blackbody with fixed yS = 2 and a best-fit 
temperature of Ti^^t = 36.5 K. We also plotted acceptable fits (see text) as dotted lines, con'espoding to Tjuji = 33 K and T^^mt = 41 K. The right panel shows 
the full IR SED including all the photome tric points in Tab l e[T](th e colour code for the dift'erent instruments appeal's in the figure inset), as well as our extracted 
IRS spectrum and the MIPS-SED data of lGallimore et al.l ilOlOl) as red solid lines. We have also plotted our modelling to the mid and fai-IR data using two 
modified blackbodies (solid black line), with fixed /? = 2 and best-fit temperatures of r^^^j = 67 K and T^'^^^ = 35 K for the warm and cold components, 
respectively (dashed lines). The grey dotted-dashed line is the fitted type 2 AGN template derived in Section lX2l 



optical spectroscopy. Table |4] summarises the relevant parameters 
resulting from the present work which, in addition to table 8 in 
ISchweizer & Seitzen ( l2007h . shows the main properties of Mrk 938. 
The analysis of the X-ray data favours the presence of a highly 
obscured AGN, where the SB emission prevails at soft energies 
and the active nucleus dominates at hard energies. The absorption- 
corrected hard X-ray luminosity is too high to be originated only in 
a pure starburst. 



We have performed a decomposition of the 5pifzer/IRS spec- 
trum and recovered the AGN contribution to the IR luminosity, 
which is overwhelmed by the intense emission of star formation. 
The best fit model to the IRS data was achieved with a combination 
of a type 2 AGN plus a 10'^ Lq starburst templates. Using the de- 
rived fractional contribution of the AGN to the luminosity at 24 \im 
we derived an AGN bolometric contribution to the IR luminosity of 
~2%, which is in agreement with previous estimations. This sup- 
ports the proposed scenario that intense dusty starbursts are respon- 
sible for the high IR luminosities in most local LIRGs. The analysis 
of IR images shows that the IR emission of Mrk 938 originates in 
a relatively compact region, with a size between 0.5 and 2 kpc. We 
have also demonstrated that most of the on-going SFR is produced 
in an obscured environment as expected for LIRGs. 



The wealth of multi-wavelength archival data of this source 
has been used in conjunction with our own observations to con- 
strain the SED. We have used Herschel imaging data for the first 
time to constrain the cold dust emission with unprecedented accu- 
racy. We have derived dust temperatures of 63 K and 35 K for the 
warm and cold components, respectively. Using the results of our 
SED fitting we have derived a cold dust mass of 3 x 10^, Mq and a 
very high mass ratio of the cold to the warm dust components. This 
demonstrates that emission of cold dust clearly dominates the IR 
SED. 
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Table 4. Summary of the properties of Mrk 938 derived in tiiis paper. 



AGN Properties 



Hard X-ray limiinosity 


i2-10keV 


1.4+0-^ X 10«ergs-i 
1.6+i-ixl0«ergs-' 
2.6+fl X lO'^^ergs'i 


absorption corrected 


Bolometric luminosity 


Lboi(AGN) 


from X-ray and bolometric correction 
from IRS decomposition 


AGN contribution 


Bolometric 


2+2% 


from IRS decomposition 




0.5 -2keV 


^5% 


from comparison with SFR-due X-ray luminosity 




2- lOkeV 


^ 85% 


from comparison with SFR-due X-ray luminosity 




24 ^im 


5 (±3)% 


from IRS decomposition 



SB Properties 



from MIPS 24 |.im, Kroupa IMF 
from IRAC 8 |xm (FWHM) 
from MIPS 24 \im, PACS 70 nm (FWHM) 
from dust mass, assuming TioOyum = 0.5-1 

two temperature fit, fixed j8 = 2 



Star Formation Rate SFRobscured 42M0yr~ 

Size IR emitting region dm. ~ 1 kpc 

2 kpc 
1 - 1.6kpc 

Dust Temperature T^^^, 35 (±4) K 

^Lt 63(±4)K 
Dust Mass Mdast 3(±l)xl0^Mo 

Dust mass Ratio ^ust^^ust '^^ 
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